Elucidating the Role of OXPHOS Variants in Asthenozoospermia: Insights from Whole Genome Sequencing and an In Silico Analysis.

Infertility is a global health challenge that affects an estimated 72.4 million people worldwide. Between 30 and 50% of these cases involve male factors, showcasing the complex nature of male infertility, which can be attributed to both environmental and genetic determinants. Asthenozoospermia, a condition characterized by reduced sperm motility, stands out as a significant contributor to male infertility. This study explores the involvement of the mitochondrial oxidative phosphorylation (OXPHOS) system, crucial for ATP production and sperm motility, in asthenozoospermia. Through whole-genome sequencing and in silico analysis, our aim was to identify and characterize OXPHOS gene variants specific to individuals with asthenozoospermia. Our analysis identified 680,099 unique variants, with 309 located within OXPHOS genes. Nine of these variants were prioritized due to their significant implications, such as potential associations with diseases, effects on gene expression, protein function, etc. Interestingly, none of these variants had been previously associated with male infertility, opening up new avenues for research. Thus, through our comprehensive approach, we provide valuable insights into the genetic factors that influence sperm motility, laying the foundation for future research in the field of male infertility.

PDE4 Gene Family Variants Are Associated with Response to Apremilast Treatment in Psoriasis.

Moderate-to-severe psoriasis (Ps) treatment includes systemic drugs and biological agents. Apremilast, a small molecule primarily metabolized by cytochrome CYP3A4, modulates the immune system by specifically inhibiting phosphodiesterase type 4 (PDE4) isoforms and is currently used for the treatment of Ps and psoriatic arthritis (PsA). Clinical trials and real-world data showed variable efficacy in response among Ps patients underlying the need for personalized therapy. This study implements a candidate-gene and a network-based approach to identify genetic markers associated with apremilast response in forty-nine Greek Ps patients. Our data revealed an association of sixty-four SNPs within or near PDE4 and CYP3A4 genes, four SNPs in ncRNAs ANRIL, LINC00941 and miR4706, which influence the abundance or function of PDE4s, and thirty-three SNPs within fourteen genes whose protein products either interact directly with PDE4 proteins or constitute components of the cAMP signaling pathway which is modulated by PDE4s. Notably, fifty-six of the aforementioned SNPs constitute eQTLs for the respective genes in relevant to psoriasis tissues/cells implying that these variants could be causal. Our analysis provides a number of novel genetic variants that, upon validation in larger cohorts, could be utilized as predictive markers regarding the response of Ps patients to apremilast treatment.

MassArray Genotyping as a Selection Tool for Extending the Shelf-Life of Fresh Gilthead Sea Bream and European Seabass.

In modern aquaculture, genomics-driven breeding programs have emerged as powerful tools for optimizing fish quality. This study focused on two emblematic Mediterranean fish species, the European seabass (Dicentrarchus labrax) and the gilthead sea bream (Sparus aurata), with a primary aim of exploring the genetic basis of white muscle/fillet degradation in fresh fish following harvest. We identified 57 and 44 missense SNPs in gilthead sea bream and European seabass, respectively, located within genes encoding for endogenous proteases responsible for fillet quality. These SNPs were cherry-picked based on their strategic location within the catalytic/regulatory domains of endogenous proteases that are expressed in the white muscle. Using MassArray technology, we successfully associated differentiated enzymatic activity of those endogenous proteases post-harvest as a phenotypic trait with genetic polymorphism of six SNPs in gilthead sea bream and nine in European seabass. These findings can be valuable attributes in selective breeding programs toward the extension of freshness and shelf life of these species. The integration of MassArray technology into breeding programs offers a cost-effective strategy for harnessing the potential of these genetic variants to enhance the overall quality of the final product. Recognizing that fresh fish perishability is a challenge, extending shelf-life is pivotal in reducing losses and production costs.

Heat Stress: A Serious Disruptor of the Reproductive Physiology of Dairy Cows.

Global warming is a significant threat to the sustainability and profitability of the dairy sector, not only in tropical or subtropical regions but also in temperate zones where extreme summer temperatures have become a new and challenging reality. Prolonged exposure of dairy cows to high temperatures compromises animal welfare, increases morbidity, and suppresses fertility, resulting in devastating economic losses for farmers. To counteract the deleterious effects of heat stress, cattl e employ various adaptive thermoregulatory mechanisms including molecular, endocrine, physiological, and behavioral responses. These adaptations involve the immediate secretion of heat shock proteins and cortisol, followed by a complex network of disrupted secretion of metabolic and reproductive hormones such as prolactin, ghrelin, ovarian steroid, and pituitary gonadotrophins. While the strategic heat stress mitigation measures can restore milk production through modifications of the microclimate and nutritional interventions, the summer fertility records remain at low levels compared to those of the thermoneutral periods of the year. This is because sustainment of high fertility is a multifaceted process that requires appropriate energy balance, undisrupted mode of various hormones secretion to sustain the maturation and fertilizing competence of the oocyte, the normal development of the early embryo and unhampered maternal-embryo crosstalk. In this review, we summarize the major molecular and endocrine responses to elevated temperatures in dairy cows, as well as the impacts on maturing oocytes and early embryos, and discuss the consequences that heat stress brings about in dairy cattle fertility.

Unveiling the Genetic Complexity of Teratozoospermia: Integrated Genomic Analysis Reveals Novel Insights into lncRNAs' Role in Male Infertility.

Male infertility is a global health issue, affecting over 20 million men worldwide. Genetic factors are crucial in various male infertility forms, including teratozoospermia. Nonetheless, the genetic causes of male infertility remain largely unexplored. In this study, we employed whole-genome sequencing and RNA expression analysis to detect differentially expressed (DE) long-noncoding RNAs (lncRNAs) in teratozoospermia, along with mutations that are exclusive to teratozoospermic individuals within these DE lncRNAs regions. Bioinformatic tools were used to assess variants' impact on lncRNA structure, function, and lncRNA-miRNA interactions. Our analysis identified 1166 unique mutations in teratozoospermic men within DE lncRNAs, distinguishing them from normozoospermic men. Among these, 64 variants in 23 lncRNAs showed potential regulatory roles, 7 variants affected 4 lncRNA structures, while 37 variants in 17 lncRNAs caused miRNA target loss or gain. Pathway Enrichment and Gene Ontology analyses of the genes targeted by the affected miRNAs revealed dysregulated pathways in teratozoospermia and a link between male infertility and cancer. This study lists novel variants and lncRNAs associated for the first time with teratozoospermia. These findings pave the way for future studies aiming to enhance diagnosis and therapy in the field of male infertility.

Recovery of Staphylococci from Teatcups in Milking Parlours in Goat Herds in Greece: Prevalence, Identification, Biofilm Formation, Patterns of Antibiotic Susceptibility, Predictors for Isolation.

The objectives of this work are (a) to describe staphylococci on the teatcups of milking parlours in goat farms and identify predictors for the presence of staphylococcal isolates on the teatcups, (b) to evaluate relationships with total bacterial counts and somatic cell counts in bulk-tank milk, and (c) to establish patterns of susceptibility to antibiotics for the staphylococcal isolates and identify predictors for the recovery of resistant isolates. In a cross-sectional study of 66 goat farms across Greece, swab samples were collected from 303 teatcups (upper and lower part) for staphylococcal recovery, identification, and assessment of biofilm formation. Details regarding health management on the farms (including conditions in the milking parlour) and the socio-demographic characteristics of farmers were collected by means of a structured questionnaire. A total of 87 contaminated teatcups (28.7%) were found on 35 goat farms (53.0%). Staphylococci were more frequently recovered from the upper than the lower part of teatcups: 73 versus 43 teatcups, respectively. After identification, 67 staphylococcal isolates (i.e., excluding similar isolates) were recovered from the teatcups; Staphylococcus aureus, Staphylococcus capitis, and Staphylococcus equorum predominated. Of these isolates, 82.1% were biofilm-forming. In multivariable analysis, the annual incidence of clinical mastitis in the herd emerged as the only significant factor associated with the isolation of staphylococci from the teatcups. Of the 67 isolates, 23 (34.3%) were resistant to at least one antibiotic, and 14 (22.4%) were multi-resistant. Resistance was found most commonly against penicillin and ampicillin (22.4% of isolates), fosfomycin (17.9%), clindamycin (14.9%), erythromycin, and tetracycline (13.4%). In multivariable analysis, the annual incidence of clinical mastitis in the herd and the use of detergent for parlour cleaning emerged as significant factors associated with the isolation of staphylococci resistant to antibiotics.

Gene-by-Sex Interactions: Genome-Wide Association Study Reveals Five SNPs Associated with Obesity and Overweight in a Male Population.

Obesity is a chronic health problem associated with severe complications and with an increasing prevalence in the Western world. Body-fat composition and distribution are closely associated with obesity, but the human body's composition is a sexually dimorphic trait, as differences between the two sexes are evident even from fetal life. The effect of sex hormones contributes to this phenomenon. However, studies investigating gene-by-sex interactions for obesity are limited. Therefore, the aim of the present study was to identify single-nucleotide polymorphisms (SNPs) associated with obesity and overweight in a male population. A genome-wide association study (GWAS) that included 104 control, 125 overweight, and 61 obese subjects revealed four SNPs associated with overweight (rs7818910, rs7863750, rs1554116, and rs7500401) and one SNP (rs114252547) associated with obesity in males. An in silico functional annotation was subsequently used to further investigate their role. Most of the SNPs were found in genes regulating energy metabolism and homeostasis, and some of them were expression quantitative trait loci (eQTL). These findings contribute to the understanding of the molecular mechanisms underlying obesity-related traits, especially in males, and pave the road for future research toward the improvement of the diagnosis and therapy of obese individuals.

Genetic diversity and thermotolerance in Holstein cows: Pathway analysis and marker development using whole-genome sequencing.

Heat stress causes extensive losses in the dairy sector, due to negative effects on milk production and reproduction. Cows have evolved a series of protective mechanisms, (physiological, biochemical, behavioural) to cope with the thermostressing environments, which have allowed the preservation of productive and reproductive potential of specific animals during summer; these animals are considered thermotolerant and could be used to design programs of selective breeding. These programs, targeting the generations of a population of heat-resistant animals, would increase the frequency of the desired phenotypes, tackling the financial losses on one hand and reducing the carbon footprints of the dairy sector on the other. The development of genomics techniques has enabled genome wide variant calling, to detect SNPs associated with the desired phenotypes. In this study, we used a comparative genomics approach to detect genetic variation associated with thermotolerance and to design molecular markers for characterizing the animals as tolerant/sensitive. A total of 40 cows from each group were split in four sequencing pools and a whole-genome sequencing approach was used. Results and conclusion: Genome-wide genetic variation between groups was characterized and enrichment analysis revealed specific pathways which participate in the adaptive mechanisms of thermotolerance, implicated into systemic and cellular responses, including the immune system functionality, Heat Stress and Unfolded Protein Response. The markers made a promising set of results, as specific SNPs in five genes encoding for Heat Shock Proteins were significantly associated with thermotolerance.

Investigating the Impact of a Curse: Diseases, Population Isolation, Evolution and the Mother's Curse.

The mitochondrion was characterized for years as the energy factory of the cell, but now its role in many more cellular processes is recognized. The mitochondrion and mitochondrial DNA (mtDNA) also possess a set of distinct properties, including maternal inheritance, that creates the Mother's Curse phenomenon. As mtDNA is inherited from females to all offspring, mutations that are harmful to males tend to accumulate more easily. The Mother's Curse is associated with various diseases, and has a significant effect on males, in many cases even affecting their reproductive ability. Sometimes, it even leads to reproductive isolation, as in crosses between different populations, the mitochondrial genome cannot cooperate effectively with the nuclear one resulting in a mito-nuclear incompatibility and reduce the fitness of the hybrids. This phenomenon is observed both in the laboratory and in natural populations, and have the potential to influence their evolution and speciation. Therefore, it turns out that the study of mitochondria is an exciting field that finds many applications, including pest control, and it can shed light on the molecular mechanism of several diseases, improving successful diagnosis and therapeutics. Finally, mito-nuclear co-adaptation, paternal leakage, and kin selection are some mechanisms that can mitigate the impact of the Mother's Curse.

Whole-Genome Profile of Greek Patients with Teratozοοspermia: Identification of Candidate Variants and Genes.

Male infertility is a global health problem that affects a large number of couples worldwide. It can be categorized into specific subtypes, including teratozoospermia. The present study aimed to identify new variants associated with teratozoospermia in the Greek population and to explore the role of genes on which these were identified. For this reason, whole-genome sequencing (WGS) was performed on normozoospermic and teratozoospermic individuals, and after selecting only variants found in teratozoospermic men, these were further prioritized using a wide range of tools, functional and predictive algorithms, etc. An average of 600,000 variants were identified, and of them, 61 were characterized as high impact and 153 as moderate impact. Many of these are mapped in genes previously associated with male infertility, yet others are related for the first time to teratozoospermia. Furthermore, pathway enrichment analysis and Gene ontology (GO) analyses revealed the important role of the extracellular matrix in teratozoospermia. Therefore, the present study confirms the contribution of genes studied in the past to male infertility and sheds light on new molecular mechanisms by providing a list of variants and candidate genes associated with teratozoospermia in the Greek population.

Oviductal epithelial cells transcriptome and extracellular vesicles characterization during thermoneutral and heat stress conditions in dairy cows.

In this study, the transcriptome of oviductal epithelial cells and certain characteristics of their extracellular vesicles of dairy cows were described under thermoneutral and heat stress conditions. Twenty cows were compared in springtime at THI = 65.6 ± 0.90 and in summertime at THI = 78.36 ± 2.73. During each season, the estrous cycles of the cows were synchronized, and on day 3 of the ensuing cycle, a blood sample was collected for progesterone determination, while their oviducts were collected after slaughter. Epithelial cells and oviductal fluid were collected from the oviduct ipsilateral and contralateral to the corpus, respectively. For the gene expression study, a comparative transcriptomic approach, using RNASeq, was performed on cells collected from the ipsilateral and the contralateral oviducts. The size and the concentration of extracellular vesicles (EVs) at both seasons were analyzed using Transmission Electron Microscopy and Nanoparticle tracking analysis and specific proteins were detected by Western blotting. Progesterone concentration was higher during the thermoneutral period. Between seasons, divergent expression of genes related to immune system, contractility, gamete protection and lncRNAs was found. The size and the concentration of the EVs did not differ between seasons, however, the concentration in the ipsilateral oviduct tended to be lower (p = 0.09) from the contralateral one in the summer, but not in the spring. Our results show for the first time that HS could be involved with alterations in the oviductal cells' gene expression and in the changes in concentration of EVs in the oviductal lumen. Our results imply that the altered oviductal environment during HS could be associated with the suppressed summer fertility in dairy cows.

The Effects of Heat Shock Protein 70 Addition in the Culture Medium on the Development and Quality of In Vitro Produced Heat Shocked Bovine Embryos.

The aims of the present study were to examine the effects of HSP70 addition in the in vitro culture medium of day 3 embryos on their developmental competence and quality. Bovine oocytes (n = 1442) were in vitro matured, inseminated and cultured for the first two days according to standardized methods. The presumptive zygotes were randomly allocated in three experimental groups: Control, C (embryos cultured at 39 °C throughout the culture period), group C41 (temperature was raised to 41 °C from the 48th to 72nd h post insemination (p.i.) and then it returned at 39 °C for the remaining culture period), and group H41 (the temperature modification was the same as in C41 and during heat exposure, HSP70 was added in the culture medium). Cleavage and embryo yield were assessed 48 h p.i. and on days 7, 8, 9, respectively and gene expression in day 7 blastocysts was assessed by RT-PCR. Blastocyst yield was the highest in group C39; and higher in group H41 compared to group C41. From the gene expression analyses, altered expression of 11 genes was detected among groups. The analysis of the orchestrated patterns of gene expression differed between groups. The results of this study confirm the devastating effects of heat stress on embryo development and provide evidence that HSP70 addition at the critical stages can partly counterbalance, without neutralizing, the negative effects of the heat insult on embryos, acting mainly through mechanisms related to energy deployment.

Heat Shock Protein 70 Improves In Vitro Embryo Yield and Quality from Heat Stressed Bovine Oocytes.

Heat shock protein 70 (HSP70) is a chaperon that stabilizes unfolded or partially folded proteins, preventing inappropriate inter- and intramolecular interactions. Here, we examined the developmental competence of in vitro matured oocytes exposed to heat stress with or without HSP70. Bovine oocytes were matured for 24 h at 39 °C without (group C39) or with HSP70 (group H39) and at 41 °C for the first 6 h, followed by 16 h at 39 °C with (group H41) or without HSP70 (group C41). After insemination, zygotes were cultured for 9 days at 39 °C. Cleavage and embryo yield were assessed 48 h post insemination and on days 7, 8, 9, respectively. Gene expression was assessed by RT-PCR in oocytes, cumulus cells and blastocysts. In C41, blastocysts formation rate was lower than in C39 and on day 9 it was lower than in H41. In oocytes, HSP70 enhanced the expression of three HSP genes regardless of incubation temperature. HSP70 at 39 °C led to tight coordination of gene expression in oocytes and blastocysts, but not in cumulus cells. Our results imply that HSP70, by preventing apoptosis, supporting signal transduction, and increasing antioxidant protection of the embryo, protects heat stressed maturing bovine oocyte and restores its developmental competence.

Developmental competence of heat stressed oocytes from Holstein and Limousine cows matured in vitro.

The negative effects of heat stress on dairy cattle's fertility have been extensively studied, but the relevant knowledge for beef cattle is rather limited. The aims of this study were to investigate the effects of HS during in vitro maturation on the developmental potential of oocytes derived from Limousine and Holstein cows and to estimate the effect of the differential gene expression of important genes in oocytes, cumulus cells and blastocysts in the growth competence between the breeds. In seven replicates, cumulus oocyte complexes from Holstein and Limousine cows were matured for 24 hr at 39°C (controls C; Hol_39, Lim_39) or at 41°C from hour 2 to hour 8 of IVM (treated T; Hol_41, Lim_41), fertilized, and presumptive zygotes were cultured for 9 days at 39°C. Cleavage and embryo formation rates were evaluated 48 hr post-insemination and on days 7, 8 and 9, respectively. From all groups, subsets of cumulus cells, oocytes and blastocysts were analysed for the relative expression of genes related to metabolism, stress, apoptosis and placentation. No difference was detected in cleavage rate or in blastocyst formation rate among the control groups. In both breeds, heat stress reduced blastocyst yield, but at all days the suppression was higher in Limousines. In Holsteins, altered gene expression was detected in cumulus cells (G6PD, GLUT1) and blastocysts (PLAC8), while in Limousines, differences were found in oocytes (G6PD, HSP90AA1), in cumulus cells (CPT1B, HSP90AA1, SOD2) and blastocysts (DNMT, HSP90AA1, SOD2). It appears that Holstein COCs are more tolerant than Limousine COCs, possibly due to compulsory, production driven selection.

HLA class II peptide-binding-region analysis reveals funneling of polymorphism in action.

BACKGROUND: HLA-class II proteins hold important roles in key physiological processes. The purpose of this study was to compile all class II alleles reported in human population and investigate patterns in pocket variants and their combinations, focusing on the peptide-binding region (PBR). METHODS: For this purpose, all protein sequences of DPA1, DQA1, DPB1, DQB1 and DRB1 were selected and filtered, in order to have full PBR sequences. Proportional representation was used for pocket variants while population data were also used. RESULTS: All pocket variants and PBR sequences were retrieved and analyzed based on the preference of amino acids and their properties in all pocket positions. The observed number of pocket variants combinations was much lower than the possible inferred, suggesting that PBR formation is under strict funneling. Also, although class II proteins are very polymorphic, in the majority of the reported alleles in all populations, a significantly less polymorphic pocket core was found. CONCLUSIONS: Pocket variability of five HLA class II proteins was studied revealing favorable properties of each protein. The actual PBR sequences of HLA class II proteins appear to be governed by restrictions that lead to the establishment of only a fraction of the possible combinations and the polymorphism recorded is the result of intense funneling based on function.

Short term temperature elevation during IVM affects embryo yield and alters gene expression pattern in oocytes, cumulus cells and blastocysts in cattle.

Heat stress causes subfertility in cattle by inducing alterations in steroidogenic capacity, follicular function and ovulation defects, which eventually negatively affect oocyte quality and embryo survival. Here, the effects of short, moderate temperature elevation during IVM, on embryo yield, and on the expression of various genes was evaluated. In 8 replicates, cumulus oocyte complexes (COCs) were matured for 24 h at 39 °C (controls n = 605) or at 41 °C from hour 2 to hour 8 of IVM (treated, n = 912), fertilized, and presumptive zygotes were cultured for 9 days at 39 °C. Cleavage and embryo formation rates were evaluated 48 h post insemination and on days 7, 8, 9 respectively. Cumulus cells, oocytes and blastocysts from 5 replicates were snap frozen for the relative expression analysis of genes related to metabolism, thermal and oxidative stress response, apoptosis, and placentation. In treated group, cleavage and embryo formation rates were statistically significantly lower compared with the control (cleavage 86.7% vs 74.2%; blastocysts: day 7, 29.9% vs 19.7%, day 8, 34.2% vs 22.9% and day 9 35.9% vs 24.5%). Relative mRNA abundance of three genes in cumulus cells (HSP90AA1, CPT1B, G6PD) and three genes in blastocysts (DNMT3A, PLAC8, GPX1) indicated significantly different expression between groups (p < 0.05)., The expression of G6PD, SOD2, GXP1 in oocytes and PTGS2 in blastocysts tended to differ among groups (0.05

Islands and hybrid zones: combining the knowledge from "Natural Laboratories" to explain phylogeographic patterns of the European brown hare.

BACKGROUND: The aim of the study was to use hybrid populations as well as island populations of the European brown hare (Lepus europaeus) to explore the effect of evolutionary events, such as the post-deglaciation translocations, spontaneous and human-mediated, local adaptation and the genetic drift in the shaping of the phylogeographic patterns of the species. For this purpose, we used molecular markers, both nuclear and mitochondrial, that are indicative for local adaptation as well as neutral markers to elucidate the patterns of population differentiation based on geographic isolation and the clade of origin. To broaden our analysis, we included data from our previous studies concerning mainland populations, to explore the genetic differentiation in the base of the geographic origin (mainland/island) of the populations. RESULTS: Our results suggest that local adaptation shapes the differentiation in both genomes, favoring specific alleles in nuclear genes (e.g. DQA) or haplotypes in mtDNA (e.g. Control Region, CR). mtDNA variation was found to be in a higher level and was able to give a phylogeographic signal for the populations. Furthermore, the degree of variation was influenced not only by the geographic origin, but also by the clade of origin, since specific island populations of Anatolian origin showed a greater degree of variation compared to specific mainland populations of the European clade. Concerning the hybrid population, we confirmed the existence of both clades in the territory and we provided a possible explanation for the lack of introgression between the clades. CONCLUSION: Our results indicate that the Quaternary's climatic oscillations played a major role in the shaping of the phylogeographic patterns of the species, by isolating populations in the distinct refugia, where they adapted and differentiate in allopatry, leading to genome incompatibilities observed nowadays.

Mitogenomic analysis in European brown hare (Lepus europaeus) proposes genetic and functional differentiation between the distinct lineages.

European brown hare is a small game species spreading across Europe to Asia Minor, with important economic traits. Population genetics studies using mitochondrial DNA markers have revealed the existence of two major phylogeographic lineages, the European and the Anatolian. European lineage is further divided in the European type halpogroup and south-eastern European type halpogroup, while Anatolian consists only by the Anatolian/Middle Eastern type halpogroup. All three haplogroups show a discrete geographical distribution, with an overlapping zone forming in North-East Greece and Bulgaria, forming a contact zone. Despite the existence of a contact zone, European haplotype was never detected in Anatolia and vice versa, proposing the presence of genetic barriers responsible for this phenomenon. In this study, we analyzed the whole mitochondrial genomes of specimens originating from both lineages, aiming to detect the genetic and functional differentiation of the oxidative phosphorylation complexes that are encoded by mtDNA that could lead gradually to the reproductive isolation of the lineages.

SNP Identification through Transcriptome Analysis of the European Brown Hare (Lepus europaeus): Cellular Energetics and Mother's Curse.

The European brown hare (Lepus europaeus, Pallas 1778) is an important small game species in Europe. Due to its size and position in the food chain, as well as its life history, phenotypic variation and the relatively recent speciation events, brown hare plays an important role in the structure of various ecosystems and has emerged as an important species for population management and evolutionary studies. In order to identify informative SNPs for such studies, heart and liver tissues of three samples from the European lineage and a three-sample pool from the Anatolian lineage were subjected to RNA-Sequencing analysis. This effort resulted in 9496 well-assembled protein-coding sequences with close homology to human. After applying very stringent filtering criteria, 66185 polymorphic sites were identified in 7665 genes/cds and 2050 of those polymorphic sites are potentially capable of distinguishing the European from the Anatolian lineage. From these distinguishing mutations we focused on those in genes that are involved in cellular energy production, namely the glycolysis, Krebs cycle and the OXPHOS machinery. A selected set of SNPs was also validated by Sanger sequencing. By simulating the three European individuals as one pool, no substantial informative-SNP identification was lost, making it a cost-efficient approach. To our knowledge this is the first attempt to correlate the differentiation in both nuclear and mitochondrial genome between the two different lineages of L. europaeus with the observed spatial partitioning of the lineages of the species, proposing a possible mechanism that is maintaining the reproductive isolation of the lineages.